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The randomly driven Navier-Stokes equation without pressure in d-dimensional space is considered 
as a model of strong turbulence in a compressible fluid. We derive a closed equation for the velocity- 
gradient probability density function. We find the asymptotics of this function for the case of the 
gradient velocity field (Burgers turbulence), and provide a numerical solution for the two-dimensional 
case. Application of these results to the velocity-difference probability density function is discussed. 



o\ 

On 
(N 

> 
O 
O 

in 
o 

00 
On 



I 

Oh: 
^ ■ 



X 



PACS Number(s): 47.27.Gs, 03.40. Kf, 52.35. Ra. 
I. INTRODUCTION 

The Burgers equation with a random external force is 
considered to be the first exactly solvable model of Id tur- 
bulence and has been extensively studied in recent years 
^]-|rT|. Though rather simplified, this model can serve 
as a test model for some general ideas within the theory 
of strong turbulence. In 1995, methods of quantum field 
theory were applied to this problem by A. Polyakov |^] 
which enabled a qualitative explanation of velocity- 
dilference probability density functions (PDFs) measured 
numerically by A. Chekhlov and V. Yakhot (|. In ||l it 
was shown that the approach [Q| allows one to obtain 
quantitatively correct results. Extensive numerical sim- 
ulations published recently by T. Gotoh and R. Kraich- 
nan show that the predictions of arc quite ac- 
curate, and coincide with the numerical simulations to 
within about 5%. V. Yakhot has shown in that the 
ideas introduced in can have much wider application, 
and can also work for incompressible velocity fluctua- 
tions. 

We believe that the operator product expansion 
(OPE), introduced in ^ to take into account the vis- 
cous term, is an adequate language to treat compress- 
ible turbulence in higher dimensions as well, where shock 
structures and associated local dissipation persist. In 
the present paper we find a closed equation for the 
velocity-gradient probability-density function (PDF) for 
compressible turbulence in any number of dimensions. 
We investigate the asymptotics of the PDF and present 
the numerical solution for the 2d case. 

The basic equation we will study is the following: 



Uf + (u • V)u = vAu + f 



(1) 



The force / is chosen to be Gaussian with zero mean 
and white in time variance: 



(r(x,i)/'(x',0> ='5(t-i')«''(x-x') 



(2) 



where the k function is concentrated at some large 
scale L, and can be expanded as follows 



for y <^ L. We assume that the steady states for velocity 
gradient and velocity difference exist; for this we can re- 
quire, for example, that periodic boundary conditions on 
a scale much larger that L be imposed, and that the zero 
harmonic in the k function be absent. These assumptions 
are usually used in numerical simulations |p|,|lo|. 

In this paper we appeal to the results obtained 
for the Id Burgers turbulence without pressure 
in In particular, we are interested in the 

velocity-gradient PDF P{du'^ / dx^) and the velocity- 
dilference PDF P-u (u(xi) — u(x2)), where the velocities 
are taken at the same time at some fixed points xi 
and X2. The physical picture presented in these papers 
allows us to consider a general phenomenon such as inter- 
mittency on a rigorous basis; it is related to the sponta- 
neous breakdown of the Galilean invariance of the forced 
equation and to the algebraic decay of the PDFs. We will 
not repeat these arguments here; instead, we will concen- 
trate on the main ideas which allow us to consider the 
multi-dimensional case. 

We will be interested in the case of small dissipation u, 
and will consider distances \xi — X2 \ <^ L. The following 
order of the limits should be considered to get the steady 
state: we first set i — > oo, and then consider the limit 



II. VELOCITY-GRADIENT PDF 

To proceed quantitatively, consider the following char- 
acteristic function [Z function) for the velocity gradient 
ul = du^dx'': 



Z{aki) = {e-Kp{iakiuf)) 



(4) 



Due to Eq. (|l|), this function satisfies the following 
Fokker-Planck equation: 



Z = i<T, 



'■^ diJikdiJii 



dZ 



(3) 



^-^CTik {cTik — CTki)] Z + D, 



(5) 
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where summation over repeated indices is assumed. For 
simplicity, we set ki = 1. To derive this equation we dif- 
ferentiated (^) with respect to t and expressed it^ using 
Eq. (|l|). We made use of the following identities: 



92 



daikddii 



- {u\ exp(icr„„M™)) = iSiko§-^Z , 

and 



(6) 



(7) 



(8) 



The D term describes the contribution of the dissipation 
and in steady state is given by: 



D ~ lim (i^icr^s exp(i(Tfc;ii; )) 



(9) 



Without this term the steady state in Eq. (|^) does 
not exist. This term can not be closed without further 
assumptions. In [Q, assumptions about scaling invari- 
ance, Galilean invariance, and the operator product ex- 
pansion were applied to close the analogous term for the 
velocity- difference PDF. It is not obvious a priori that 
these methods can be applied to our problem, since the 
limits y — > and 1/ — > may not be interchangeable. 

Nevertheless, it was observed in the numerical sim- 
ulations in that the velocity-difference and the 
velocity-gradient PDFs coincide for the Galilean invari- 
ant region Au <C Urms (urms — (KoLy/^) in the 
one-dimensional case. This suggests that the velocity- 
difference PDF is contributed to by smooth parts of the 
velocity field in the Galilean invariant region, and there- 
fore the limits y and v ~* are interchangeable. 

Another important result of is that the (3 anomaly 
introduced in [Q is absent for regular forcing. We as- 
sume that this is true for the multi-dimensional case as 
well. Under this assumption, the D term in the multi- 
dimensional case should be expanded as 



D = aZ 



(10) 



This is the only assumption we use in what follows. We 
refer the reader to papers ||l|-||,^ for more details and 
discussions on the underlying ideas. We will see that this 
assumption is self-consistent; the anomaly a can be found 
from the conditions of positivity, finiteness, and normal- 
izability of the PDF. These conditions can be easily im- 
posed on the PDF in u-space. We therefore transform 
Eq. (0) to M-space, using: 



P{Vu) = / daZ{a)e 



where dcr = ^ dcTifc is the measure in d 
space of the elements of the matrix Uik- 
state, the equation takes the form: 



(11) 

^-dimensional 
In the steady 



0^ 1 l+g d ( _d_ I d 



+ 



+ 



i-a d ( _a a_ 

2 du^^ \ du^^ du^ 



P 



-aP. 



(12) 



This is the general equation for the PDF. The force in this 
equation can have different symmetry properties, which 
correspond to different values of the parameter a in (|^) . 
No restrictions have so far been imposed on the velocity 
field either. 

Eq. ( I12I) can be simplified for the gradient force f — 
\7(j). This choice corresponds to a = 1 and allows us to 
look for a solution in the factorized form: 



P 



.i<k 



P{u) , 



(13) 



where ii denotes the symmetric part of the matrix Uik- 
Physically, this means that we have restricted our con- 
sideration to gradient fluctuations of the velocity field, 
u = V/i. We refer to this case as multi-dimensional Burg- 
ers turbulence. It was considered by completely different 
methods in [^,0. The spirit of our method is most close 
to the consideration of Eq. (12) with ansatz ( p^ ) can 
be cast into the following form: 



(d- 



2)u\P + ulu^^P- 



p 



aP. 



(14) 



In what follows we will consider only the function P, and 
will omit the overtilde sign. 

Eqs. ( p^ and (jl^) help to reveal the physical sense of 
the a anomaly. Integrating these equations with respect 
to ul one gets: 



-a, 



(15) 



which means that this anomaly describes the average 
measure loss due to compressibility and presence of 
shocks. It can also be interpreted as the mean rate of 
density accumulation on shocks in the Lagrangian pic- 
ture: 



p{y,t) + ul p{y,t) + a p{y,t) 



0, 



(16) 



where y is the Lagrangian coordinate, u\ represents the 
smooth part of the velocity field, and a = (ul) shocks- 
The Lagrangian PDF P{u%,p;y) and the Eulerian 
PDF P(u^;x) are related as follows: 



Piul^x) 



J ^Piul,p;y), 



(17) 
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that leads to Eq. (|12D. We note an interesting analogy 
with a mean field approximation: a is introduced as the 
mean field in the dynamical equations and then is found 
self-consistently firom Eq. ( [l^ ) . The interpretation ( [l6| ) is 
important, since it allows one to introduce the anomaly 
on the level of the stochastic Langevin equation. 

In general, the PDF should depend only on invariants 
with respect to space rotations. For the d-dimensional 
space, there are exactly d such invariants, which can be 
chosen as the eigenvalues of the matrix Uik- Let us de- 
note them as Ai, A2, . . . , \d- Eq- (0) can be rewritten 
for the function P depending on only these variables: 



E {d + 2)Xk + X 

k=l 
d 



2 d 
kd\k 



P 



a 



k=l 
_ _d_ 



d 



P 



P 



-aP. (18) 



To derive this equation we used the following expression 
for the matrix Laplacian, known in the theory of matrix 
models mjll: 



^ dX 



aP_ 1 1 f d d_ 

■i.k 



(19) 



Equation (|lj) has an infinite number of solutions. The 
physically reasonable solution should satisfy conditions of 
positivity, finiteness, and normalizability; exactly in the 
same manner as the ground state is determined in quan- 
tum mechanics. The solution should also be symmetrical 
with respect to the arguments Ai, . . . , A^. These condi- 
tions should determine the unknown parameter a. This 
parameter depends only on the symmetry properties of 
the external force and on the space dimensionality. 

As in the one-dimensional case, the asymptotics of the 
solutions can be found by balancing different terms in 
eq. (p^. If we balance the advective and force terms, we 
will get the PDF tail in the region where the dissipative 
effects are negligible. In the one-dimensional case this 
corresponds to the right tail of the PDF. This tail de- 
cays hyper-exponentially fast. In the multi-dimensional 
case the analogous asymptotic should have the form: 
P (X exp(S'(Ai, . . . , Arf)). The function 5* should be sym- 
metric with respect to its arguments Ai,...,Ad. The 
asymptotic can be simply found for large positive A's in 
the direction close to Ai = . . . = Xd: 



P oc exp 



-A3 



(20) 



where A = Ti{uik) = Ai -I- . . . -|- A^. The same asymptotic 
for large A can also be obtained by the instanton meth- 
ods applied directly to quantum mechanics (p^). 

The tail, corresponding to large negative A's (the "left" 
tail) decays rather slowly. The explanation is simple. 
Burgers shocks always have negative velocity jumps, and 
therefore large positive velocity gradients are less proba- 
ble than large negative ones. The left tail is determined 



by large negative gradients, and to obtain it we should 
neglect the force term in (|l^). We find: 



P (X 



G{Iu 



(Ai...Ad) 



d+2 



(21) 



where G is some function, and li^ = {Xi — Xk)/XiXk 
are d(d — l)/2 invariants of the characterictic equations 
for (|l8|). The finite solution which is non- vanishing for 
Ai = . . . = Arf has the form: 



P (X 



1 



ix,...Xdr+^ 



Det-(''+2)({t,fc) 



(22) 



The function obtained from ( p^ ) should be normal- 
ized with respect to the flat measure in the d{d + l)/2- 
dimensional space of elements of the symmetric ma- 
trix Uik- In A-space, this normalization is performed as 
follows: 



/ P(A)|A(A)indAfe = l, 

•' k=l 



(23) 



where A(A) = Yli^j {Xi — Xj) is the VanderMonde deter- 
minant; for details, see p^-p^. 



III. 



NUMERICAL SOLUTION FOR THE 
TWO-DIMENSIONAL CASE 



In this section we solve Eq. (^8|) numerically in the 
two-dimensional case. The purpose of these calculations 
is to show that Eq. (||) with the anomaly term ( p^ ) does 
have a steady state, at least for the gradient velocity field. 

We have used the relaxation method and started with 
some arbitrary but symmetrical initial distribution. The 
numerical value for the anomaly turned out to be a = 
1.30 ± 0.02. The PDF has hyper-exponential and power- 
like tails and is presented in Fig. |l|. The PDF is normal- 
ized according to (p3|). Plotted on the horizontal axes 
are Ai and A2. 




FIG. 1. Velocity-gradient probability-density func- 
tion P(Ai, A2). 
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Fig. g shows the same PDF for the diagonal direction 
Ai = A2. The left tail decays as 1/A^, the right tail 
asymptotic is P cx exp(— A'^/48), in agreement with ( |20| ) 
and (H). P(A) is plotted vs. ^Af+Af = A/V2. 




FIG. 2. Velocity-gradient PDF P(A) for the diago- 
nal direction Ai = A2. 



IV. CONCLUDING REMARKS 

The crucial assumption in our treatment of the dis- 
sipative anomaly is the assumption that only smooth 
parts of the velocity field contribute to the velocity- 
difference PDF. After the velocity-gradient Z func- 
tion is found, the velocity-difference Z function can 
be constructed as follows: 



(24) 



i.e., we simply changed aik — > Qyk in (j^- The Fourier 
transform with respect to C will then give the velocity- 
difference PDF. As an example, let us consider the lon- 
gitudinal velocity-difference PDF: 



Pi{Au,y) oc 
/ d-UifcdCP(A) exp{i(yninkUik 



(25) 



where rii is a unit vector in the direction of y^. 
Since Pi{Au,y) does not depend on n^, one can aver- 
age with respect to all possible directions of this vector 
and get the following result: 

Pi{Au,y) oc 

i / P(A)| A(A)|<5 (A.n2 - ^) ^(1 - n^) f[ dA.dn,. (26) 

In the two-dimensional case this integral can be sim- 
plified further: 



Pi{Au,y) cx 



1 



Au/y 



dAidA2|Ai-A2|P(A) 



Au/y 



(If - ^0 {'^ - ^ 



1/2 



(27) 



In general, Pi{Au,y) can be represented as Pi{Au,y) = 
w{Au/y)/y. If we assume that for large negative Au/y 
the integral (|2q) is contributed to by the force- free 
asymptotic (|2l[)7we immediately get the left tail asymp- 
totic for the longitudinal velocity-difference PDF: 



(28) 



Analogously, one can obtain a PDF for V • u. For this 
purpose one should set SikC- Such a PDF was 

investigated numerically in [ p^ , though the Re number 
was not large enough to obtain the inertial range. 

Finally, we would like to note that the absence of the 
(3 anomaly, that we assumed in our consideration, can be 
not a universal fact. It was conjectured in ||] that dif- 
ferent dissipative regularizations (e.g. hyper-dissipation 
{—lyd^P / dx^^) can lead to different steady states. This 
assumption is natural for the language of the OPE: differ- 
ent dissipative operators should have different expansion 
coefficients a and b (we use the notation of [Q). More- 
over, some analog of the /3 anomaly can also be present 
in Eq. (|2[), since it describes a general velocity field, 
without "gradient" restriction (p^. 

These questions are under consideration, the results 
will be reported elsewhere. 
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to thank T. Gotoh, V. Gurarie, and R. Kraichnan for 
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on both the physics and the numerics of the problem, 
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